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t — ■ Abstract 

0^ , In light of the first release of a high precision measurement on the positron 

fraction by the AMS-02 experiment, we perform a detailed global analysis on the 

^^ \ interpretation of the latest data of PAMELA, Fermi-LAT, and AMS-02 in terms 

of dark matter (DM) annihilation and decay, using the GALPROP code. The 
allowed regions for the DM particle mass and its annihilation cross section or 
decay life-time are obtained for channels with leptonic final states: 2e, 2/x, 2r, 4e, 
4/x and At. We show that for DM annihilation into 2/x (4//) final states, the AMS- 
02 data alone favour a DM particle mass ~ 460 (900) GeV, which is significantly 
different from that favoured by the Fermi-LAT data on the total flux of electrons 
and positrons. The allowed regions by the two experiments do not overlap at 
99.99999% C.L., which is insensitive to the changes of the normalization and slope 
of the conventional astrophysical background. Consistent fits are obtained for 
annihilation channels with 2r (At) final states, which favours TeV scale DM particle 
with large cross sections above ~ 1.7 (4.0) x 10 -23 cm 3 s . For all the considered 
leptonic channels, the current data favour DM annihilation over decay, and show 
no significant indication of charge asymmetric DM decay. 
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1 Introduction 

Compelling evidence from astronomical observations has indicated that dark matter con- 
tributes to nearly 26% of the energy density of the Universe [JJ. Popular DM candidates 
such as the weakly interacting massive particles (WIMPs) are expected to annihilate 
or decay into standard model (SM) final states in the Galactic halo and beyond, which 
may leave imprints in the fluxes of cosmic-ray particles, including electrons, positrons, 
antiprotons, and cosmic gamma-rays. The ongoing satellite-borne experiments, such 
as PAMELA, Fermi-LAT, and AMS-02 etc. are searching for such potenital indirect 
signatures of DM with high accuracy. 

Significant progresses have been made in the recent years. For instance, the PAMELA 
collaboration has reported that the ratio of the positron flux to the total flux of electrons 
and positrons rises with increasing energy in the range 10-100 GeV, which is not expected 
from conventional astrophysical backgrounds [2j[3]. This result was later confirmed by 
Fermi-LAT, which further showed that the positron fraction continues to rise in the 
higher energy range 100-200 GeV pE]. The total flux of electrons and positrons measured 
by the ballon-borne experiment ATIC and BBP-BETS also showed an excess in 300-700 
GeV with a peak located at around 600 GeV 0(6]. Although the ATIC/BBP-BETS 
"bump" was not confirmed by Fermi-LAT [TIE] and HESS [H[T0], the featureless power- 
law spectrum measured by Fermi-LAT corresponds to a power index 3.08 [UJ which is 
harder than what expected from the convention astrophysical background, and may also 
require an exotic source of cosmic-ray electrons/positrons. 

Recently, the Alpha Magnetic Spectrometer (AMS-02) collaboration has released the 
first measurement of the positron fraction based on the collected 6.8 x 10 6 events of 
electron and positron with unprecedented accuracy [12J. The result shows a steadily 
increasing of positron fraction from 10 to ~ 250 GeV, which is consistent with the 
previous measurements by PAMELA and Fermi-LAT. The spectrum measured by AMS- 
02 is slightly lower than PAMELA for electron energy larger than ~ 40 GeV, and the 
slope of the positron fraction spectrum decreases by an order of magnitude from 20 to 
~ 250 GeV. Furthermore, the current AMS-02 data show no significant fine structure or 
anisotropy in the positron flux. 

The rising spectrum of positron fraction may have astrophysical origins, such as from 
nearby pulsars [HI [14] and supernovae remnants [TBI [16] . These explanations will be 
constrained by the anisotropy and the high energy behaviour of the positron flux. Halo 
DM annihilation or decay can provide an alternative explanation. Given the precisely 
measured shape of the positron fraction spectrum by AMS-02, more insights on the 
nature of DM can be obtained. Although at this stage it is still impossible to distinguish 
different type of exotic sources, this precision measurement on the positron fraction 



spectrum may shed new light on the origin of high energy cosmic-ray positrons. 

In this work, we perform an updated global analysis on the DM interpretation of the 
current measurement of cosmic-ray electrons and positrons, for both annihilation and 
decay scenarios. We calculate propagations of cosmic-ray particles and the astrophysical 
backgrounds using the numerical package GALPROP. Typical DM annihilation and de- 
cay channels such as 2e, 2p, 2r, 4e, 4/i and 4r are included. For decay scenario, both the 
charge symmetric and asymmetric cases are considered. The results show that for DM 
annihilation into electron and muon final states, the parameter regions determined by 
AMS-02 is highly inconsistent with that favoured by Fermi-LAT data on the total flux 
of electrons and positrons, which is insensitive to the normalization and slope of the con- 
ventional astrophysical background. Consistent fits are obtained only for r-lepton final 
states, which favours TeV scale DM with large cross sections above ~ 1.7 x 10~ 23 cm 3 s _1 . 
The current data favour DM annihilation over DM decay, and do not show indications 
of charge asymmetric decay. 

2 Sources and propagation of cosmic-ray particles 

In the diffusion models of cosmic-ray propagation, the Galactic halo within which the 
diffusion processes occur is parametrized by a cylinder with half-height Z^ and radius 
R. The densities of cosmic-ray particles are vanishing at the boundary of the halo. The 
processes of energy losses, reacceleration, annihilation, as well as the secondary sources 
of cosmic-rays are confined within the Galactic disc. The diffusion equation for the 
cosmic-ray particles is given by 

3ih 3 3 3 r v 
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where ip(r,p,t) is the number density per unit of total particle momentum which is 
related to the phase space density f(r,p,t) as ip(r,p,t) = 4np 2 f(r,p,t). For steady- 
state diffusion, it is assumed that 3ip/3t = 0. The spatial diffusion coefficient D xx is 
parametrized as: D xx = (3Do(p/po) s , where p = p/(Ze) is the rigidity of the cosmic- 
ray particle, and 5 is the power spectral index which can take different values Si or # 2 
when p is below or above a reference rigidity po- The parameter Dq is a normalization 
constant, and (3 = v/c is the velocity of the cosmic-ray particle. The values of D 
and S are determined from the measurements of the fluxes of other cosmic-ray species 
such as the ratio of Boron to Carbon (B/C) and that of isotopes 10 Be/ 9 Be, etc.. The 
convection term is related to the drift of cosmic-ray particles from the Galactic disc due 



to the Galactic wind. The direction of the wind is usually assumed to be along the 
^-direction perpendicular to the galactic disc, and is a constant Vc = [29(z) — 1]V C . The 
diffusion in momentum space is described by the reacceleration parameter D pp which 
is related to the Alfven speed V a of the disturbances in the hydro dynamical plasma as: 
Dpp = 4VaP 2 /(3D xx 5(4: — 5 2 )(4 — 5)), The momentum loss rate is by ft, and r/ and r r 
are the time scales for fragmentation and radioactive decay, respectively. 

The injection spectrum of primary electrons or nuclei is assumed to have a broken 
power low behaviour dq{p)/dp ~ p~ 7 , with 7 = 71(72) for the rigidity p below (above) 
a reference rigidity p s . Secondary cosmic-ray particles are treated as the final product 
of decay of charged pions and kaons created in collisions of primary cosmic-ray particles 
with interstelar gas. 

In this work, we solve the diffusion equation and calculate the energy losses for 
electrons by ionization, Coulomb interactions, bremsstrahlung, inverse Compton, and 
synchrotron using the numerical code G ALPROP v54 [T7H2T] which makes use of realistic 
astronomical information on the distribution of interstellar gas and other data as input, 
and consider various kinds of data including primary and secondary nuclei, electrons 
and positrons, 7-rays, synchrotron and radiation etc. in a self-consistent way. Other 
approaches based on simplified assumptions on the Galactic gas distribution which allows 
for fast analytic solutions can be found in Refs. [221426] . 

The electron/positron flux is related to the density function as $ e ± = vtp(p) e ±/(4n). 
For high energy electrons/poistrons v ~ c. At the top of the atmosphere (TOA) of 
the Earth, the fluxes of cosmic-rays are affected by solar winds and the helioshperic 
magnetic field. This effect is taken into account using the force-field approximation 
|27j . The electron/positron flux at the top of the atmosphere of the Earth $ T _P A which 
is measured by the experiments is related to the interstellar flux as: $JP a (Ttoa) = 
(Ttoa/T) 2 § e ±(T) , where Ttoa = T — 0^ is the kinetic energy of electrons/positrons at 
the top of the atmosphere of the Earth. We take (pp = 0.55 GV in numerical analysis. 
As we are interested in electrons/positrons fluxes at energies above 10 GeV, the effect 
of solar modulation is less significant. 

2.1 Primary sources of positron from DM annihilation/decay 

The primary source term from the annihilation of Majorana DM particles has the form 

ge(r ' p)= 2^ M 5>^^> ( 2 ) 

where (av) is the velocity-averaged DM annihilation cross section multiplied by relative 
velocity (referred to as cross section), and p(r) is the DM energy spatial distribution 



function. dNe /dp is the injection spectrum from DM particles annihilating into e^ via 
all possible intermediate states X with r]x the corresponding branching fractions. 

In the case of DM decay, if the DM particle is not its own antiparticle, its decay into 
charged leptons can be charge asymmetric, for instance, \ ~^ e + +Y~ while x ~^ e~+Y + , 
where Y ± can be the SM gauge boson W ± or charged Higgs boson H ± or other charged 
particle in new physics models. In the generic case where the energy density of the relic 
DM particles is also asymmetric, i.e., p x (r) 7^ p x (r), the source term can be written as 

1 \ P^ hx^r dN * X) r Q \ 

Se ± (r,p) = ^(l±e)X>— , (3) 

where p(r) = p x (r) + p x (r), e = (p x (r) - p x (r))/(p x (r) + p x (r)), and r is the life- 
time of the DM particle. The case where e = 1 (—1) corresponds to the DM particle 
decaying into e + (e~) only, and e = corresponds to the decay equally into e + and 
e~ . The phenomenology of charge asymmetric decay has been investigated previously in 
Refs. [28H3Q]. 

The injection spectra dNe /dp from DM annihilation and decay are calculated using 
the numerical package PYTHIA v8.175 [31]. For decay scenario, we assume that X^ 1 are 
much lighter than the DM particle, and neglect its mass effect in the kinematics of DM 
decay. 

The fluxes of cosmic-ray electrons and positrons from DM annihilation depend only 
weakly on the DM profile. In this work, we shall take the Einasto profile [32] 
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with oie ~ 0.17 and r s ps 20 kpc. The local DM energy density is fixed at p Q = 0.43 [33] . 
The background spectra of primary electrons $ e 5 and secondary positrons $ e + are 
calculated by GALPROP v54 [T7H2T] . The configuration of the propagation parameters 
are chosen corresponding to the "conventional" configuration [21]. The main parameters 
are: R = 20 kpc, Z h = 4.0 kpc, D = 5.75 x 10 28 cmV 1 , p = A GV, <J 1(2) = 0.34(0.34), 
V a = 36 km/s, p s = 4 GV, and 71(2) = 1.6(2.5). The interstellar magnetic field used 
in the GALPROP code is B(r,z) = B exp(—(r — r Q )/ro) exp(— |z|/-Zo) with B = 5 x 
10~ 10 Tesla, r = 10 kpc, and z = 2 kpc. In solving the propagation equation, we use 
a very fine grid in kinetic energy space with a logarithmic scale factor 1.02 in the range 
10 MeV - 10 TeV. The calculations of the electron propagation in this work are cross 
checked by using the GALPROP webrun 



3 Data selection and fitting schemes 

In order to avoid uncertainties caused by solar modulation, we consider the latest preci- 
sion cosmic-ray data from satellite-borne experiments with electron energy above 10 GeV 
which includes: 7 data points of the positron fraction from PAMELA in 2010 [3], 10 data 
points of the positron fraction from Fermi-LAT data in 2011 [3], and 42 data points of 
the positron fraction from AMS-02 [12]. For the total flux of electrons and positrons we 
consider the updated data of Fermi-LAT in 2010 which have 37 data points [8]. Thus in 
total 96 data points are included in the global fits. Some previous global fits to the earlier 
data can be found in Refs. [35H37] . Note that the electron spectrum of the Fermi-LAT 
2010 data is smoother than the previous one reported in 2009 (7], which may result in 
changes in best fitted DM particle mass and annihilation cross section or decay lift-time. 
In this work, the DM particle annihilating and decay into two-body and four-body 
charged leptonic final states e, /i, and r are considered. The relevant quantities are 
determined through \ 2 fits to the data. The expression of x 2 is given by 

( fth /exp\2 

x 2 - e ( ' ;/• } ■ p> 
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where ff 1 are the theoretical predictions for the positron fraction or total flux of electrons 
and positrons, and f^ xp are the corresponding experimental data with error ctj. The index 
% runs over all the available data points. 

In comparisons between the experimental data and theoretical predictions, we take 
into account the effect of finite energy resolution of the detectors, namely, the predicted 
fluxes are convoluted according to the energy resolution for each experimental detector. 
The energy resolution of PAMELA is nearly a constant ~ 5% above 10 GeV [38]. The 
energy resolution of Ferm-LAT is ~ 6% at 7GeV and ~ 15% at 1 TeV, and we take the 
numerical values from Ref. [8] in the calculation. For the resolution of AMS-02 detector, 
we use the parametrization a(E)/E = [(0.104/ v /£'/GeV) 2 + (0.014) 2 ] 1/2 which reaches 
~ 1.4% at high energies [T2] . 

The uncertainty in the normalization of the astrophysical background is taking into 
account by a scaling factor n. The uncertainty in the slope of the background is 
parametrized as (E/Eq) s with Eq = 1 GeV a reference energy. Thus the expressions 
for the positron fraction and total flux of electron and positron are modified as 
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where ^i 4 and $ ± are the fluxes from DM annihilation/decay and background calcu- 
lated from the GALPROP code, respectively. For the fit to each final states, we first vary 
the value of k in the fit to all the data to determine the best fitted value of k and them 
fix it in calculating the contours of allowed regions of parameters for a given confidence 
level. The cases with 5 = and 5^0 are investigated separately. 

4 Results 

We first consider the DM annihilation into two charged leptons 2e, 2/x, and 2r. The 
allowed regions in the [m x , (crv)) plane at 99% CL. corresponding to A% 2 = 9.21 for 
two variables are shown in the left column of Fig. [TJ Unlike the PAMELA positron 
fraction data which can only determine the annihilation cross section as a function of 
DM particle mass, the high precision AMS-02 data actually constrain both the mass and 
annihilation cross section in a relatively small region. For 2e final states, the allowed 
values are (av) w 4.8 x 10 -25 cm 3 s _1 and m x « 309 GeV. For 2/i final states, the 
allowed cross section and DM mass are larger (av) ~ 1.9 x 10 -24 cm 3 s _1 and m x w 458 
GeV. For It final states, the allowed DM mass is around 1.4 TeV with cross section 
(av) « 1.7 x 1(T 23 cm 3 s -1 . 

As it can be seen clearly from Fig. [U given the conventional pre- Fermi background 
calculated from GALPROP [21], at 99% C.L., neither the cross section nor the DM 
particle mass determined by the AMS-02 data is consistent with Fermi-LAT data for 
electron and muon channels. For instance, for 1\x channel, the Fermi-LAT data alone 
favour m x wl.5 TeV which is far away from the value ~ 460 GeV determined by AMS- 
02. The best fitted cross section from Fermi-LAT is (av) ~ 5 x 10 -24 cm 3 s _1 which is also 
larger roughly by a factor of two. In Fig. El detailed comparisons between the AMS-02 
and Fermi-LAT favoured regions at higher confidence levels 99.999% (99.99999%) C.L. 
corresponding to A% 2 = 23.0 (32.2) are shown. One sees that even at 99.99999% C.L. 
the two experimental results of AMS-02 and Fermi-LAT cannot be reconciled. 

We find that in the above mentioned channels, the best fitted background factors are 
similar k ~ 0.83 — 0.86 for fixed 5 = 0. The favoured small k < 1.0 is driven mostly by 
the Fermi-LAT data, as at relatively low energies the total flux of electrons and positrons 
measured by Fermi-LAT is slightly below the conventional background. We also perform 
other sets of fit with both k and 5 as free parameters. The contours of allowed region in 
(ft, 5) plane for 2e, 2/j,, 4e and 4// channel are shown in Fig. [31 The best-fitted value of 5 
is found to be compatible with zero, and at 95% C.L, the allowed range is quite narrow 
—0.04 < 5 < 0.04. These results indicate that such a large discrepancy between AMS-02 
and Fermi-ALT are unlikely to be removed by varying the normalization and slope of 
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Figure 1: Allowed regions in (m x , (cry)) plane at 99% C.L. by the data of PAMELA 
[3], Fermi-LAT (4j[H] and AMS-02 [12] for the case of DM annihilation. The first row 
correspond to DM annihilation into 2e and 4e. The second and third row correspond 
to 2/i, 4/i and 2r, At respectively. The AMS-02 allowed region is also shown in each 
sub-plot. The best fitted value of background factor k is shown in each plot. 
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Figure 2: Allowed regions by the data of AMS-02 [12] and the data of Fermi-LAT for 
the total flux of electrons and positrons [8] in (m x , (crv)) plane at 99%, 99.999% and 
99.99999% C.L. for the case of DM annihilation. The upper row correspond to DM 
particle annihilating into 2e and 4e, and the lower row correspond to the annihilation 
into 2/i, 4/x, respectively. The AMS-02 allowed region is also shown in each sub-plot. 



the conventional astrophysical background. 

In the framework of DM annihilation, the AMS-02 result is consistent with Fermi- 
LAT only for 2r and 4r channels. The combined results point to a heavy DM particle 
~ 1.4 (3.1) TeV and annihilation cross section ~ 1.7 (4) x 10~ 23 cm 3 s _1 for 2r(4r) final 
states. The required boost factor Bp which is defined as the enhancement of the DM 
annihilation cross section relative to the typical WIMP cross section, i.e., Bp = (av)/(3 x 
10~ 26 cm 3 s _1 ) reaches (9(1000), which is obviously too large, and is severely constrained 
by the nonobservation of gamma-rays from various sources in the sky. For instance, the 
constraint from the Galactic diffuse gamma-ray can reach ~ 1 x 10~ 23 cm 3 s _1 at m x wl.5 
TeV for 2r channel with isothermal DM profile [32] • Similar discussions only considering 
the AMS-02 data can be found in Refs. [401141] . 

The spectrum from the four lepton channels are slight smoother, and the fit results 
prefer larger DM masses. For instance, the best-fitted DM mass is ~ 900 GeV for 4/x 
and ~ 3 TeV for 4r final states. The required cross sections are also large, as it is shown 
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Figure 3: Contours of allowed regions in (k, 6) plane, at 95% C.L., from global fits to 
the data of PAMELA [3], Fermi-LAT gJE] and AMS-02 [12] for annihilation channels 
with 2e, 2/x, 4e and 4^u final states. 



in the right column of Fig. [TJ For the combined fits including all the experimental data, 
due to the significantly higher precision of the AMS-02, the allowed regions are almost 
identical to that determined from the AMS-02 data alone. 

The best fitted parameters as well as the contribution to the total x 2 from each 
experiment are listed in Tab. [TJ The qualities of fit are not good for final states with 
electrons and muons. In the 2/j, channel, the contributions to x 2 from AMS-02 and Fermi- 
LAT are 74.8 and 160.9, respectively. One source of the inconsistency is related to the 
incompatibility between AMS-02 and Fermi-LAT. Another reason for high x 2 is that the 
spectral feature of the measured positron fraction from AMS-02 does not match precisely 
the predictions from DM annihilation. The spectrum measured by AMS-02 features a 
visible "step" at ~ 100 GeV, which is difficult to match for a smooth continuum flux 
predicted from theoretical calculations. Among all these channels, only the 2r and At 
channels have x 2 /d.o.f < 2. 

We proceed to perform global fits for the DM decay scenario. Although the injection 
spectrum from a decaying DM particle with mass m x should be the same for the annihi- 
lating DM particle with mass m x /2, the final flux after propagation is slightly different 
due to the different DM density distribution dependences, namely, the source term is 
proportional to p(r) in the case of decay, while it is proportional to p 2 (r) in the case 
of annihilation. For Einasto profile, the final electron/positron flux after propagation 
tends to be slightly shaper in the case of DM decay. Such a small difference between de- 
cay and annihilation in the spectral feature is unlikely to be distinguished by PAMELA 
and Fermi-LAT data, but can lead to significantly different results when fitting the high 
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mode 


m x (GeV) 


B F 


XpAM 


^Fcr,e+e - 


A-Fer,e+ 


Xams 


Xt 2 ot/ dof 


2e 


309(309,315) 


16(16,17) 


6.0 


238.3 


41.1 


142.9 


428.3/94 


2/i 


458(399,527) 


64(50,81) 


5.0 


160.9 


40.9 


74.8 


281.5/94 


2t 


1446(1259,1694) 


560(438,740) 


3.6 


84.3 


43.0 


41.2 


172.1/94 


4e 


361(309,407) 


17(13,22) 


6.0 


197.3 


39.9 


137.7 


380.9/94 


4/x 


899(783,1033) 


131(103,169) 


4.4 


126.3 


41.5 


51.1 


223.3/94 


4r 


3068(2619,3667) 


1320(1007,1794) 


3.4 


79.4 


43.7 


44.6 


171.1/94 



Table 1: Best-fitted values of DM particle masses, boost factors, and \ 2 values for two 
lepton and four lepton final states in the case of DM annihilation. The allowed ranges 
for m x and Bp at 99% CL. are shown in the parentheses. The contributions to the 
total x 2 from each experimental data such the positron fraction from PAMELA (Xpam)> 
Fermi-LAT (xF ere +) and AMS-02 (xams) and the total flux of electrons and positrons 
from Fermi-LAT (XF er e+e-)- are a ^ so shown. 

precision AMS-02 data. 

We first consider the case with charge symmetric decays, namely, e = 0. We find 
that in general the fits in decay scenario lead to much larger \ 2 than that in the case of 
DM annihilation, especially for the channels with electron final states, the x 2 can reach 
~ 2000, indicating a rather poor fit. This is partly due to the previously mentioned fact 
that the current AMS-02 data feature a visible "step" at ~ 100 GeV. The spectrum after 
~ 100 GeV is suddenly lower than that expected from a smooth extrapolation from low 
energy region, which makes it more difficult to match the predicted positron spectrum 
in the case of DM decay, as it has a larger slope compared with that in the case of DM 
annihilation. We thus only show the allowed regions for fits with \x and r final states in 
Fig. HJ From the figure, one sees that there is no overlapping region between Fermi-LAT 
and AMS-02 favoured region for all the final states. The best fitted parameters as well as 
the contribution to the total \ 2 from each experiment for the case of charge symmetric 
decay are listed in Tab. |2j It is clear that 2/i and 4/i channels are not favoured by the 
AMS-02 data. In all the channels, one finds x 2 /d.o.f > 2, which means that the DM 
decay scenario is disfavoured for these final states. 

In the case of charge asymmetric decay (i.e. e^O), the predicted positron fraction is 
allowed to vary according to e without changing the prediction for total flux of electrons 
and positrons. For e = 1, the predicted positron fraction can be enhanced by a factor of 
two, compared with the case in which e = 0. On the other hand, e < will suppress the 
positron fraction. This leads to more freedom to fit the data. Note that changing the 
value of e will only modify the normalization of the positron fraction spectrum, it has 
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Figure 4: Allowed regions in (m x ,T) plane at 99% C.L. by the data of PAMELA [3], 
Fermi-LAT [HE] and AMS-02 [12] for the case of symmetric DM decay, i.e., e = 0. 
The upper row correspond to DM decay into 2/i and 4/a final states, and the lower row 
correspond to 2r, 4r final states, respectively. The AMS-02 allowed region is also shown 
in each sub-plot. 

no effect on the shape of the spectrum. 

We perform a scan over the allowed values of e through fitting to the data of 
PAMELA, Fermi-LAT, and AMS-02 for DM particle decaying into 2/i, 4/x, 2r, and 
4r final states. The background factors are fixed at k = 0.85 and 5 = 0. The resulting 
X 2 curves are shown in Fig. [5j For all these four decay channels, there is no significant 
indication that a nonzero e is favoured. This result can be understood from comparing 
the allowed regions for DM particle mass m x and life-time r by AMS-02 and Fermi-LAT 
in Fig. HI Although the values of m x favoured by the two experiments are significantly 
different, the favoured values of r are quite similar for all the decay channels. A non- 
vanishing e cannot change the best-fitted values of m x as they are determined by the 
spectral feature of the observed electron/positron fluxes, but will enlarge the difference 
in the favoured value of life-time r between AMS-02 and Fermi-LAT, which leads to the 
increase of \ 2 an d is disfavoured. 
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mode m x (GeV) r(xlO ifa s) Xpam xL,e+e- xL,e+ Xams Xtot/dof 

2/i 415[383,449] 6.6[6.2,7.0] 

2r 1209[1096,1336] 2.3[2. 1,2.5] 

Afi 814(767,881] 6.1[5. 7,6.5] 

At 2467[2235,2779] 2.1 [1.9,2.2] 



Table 2: Best-fitted values of DM particle masses, life-times, and \ 2 values for DM 
particle decaying into /i- and r-lepton final states. The allowed ranges for m x and r at 
99% CL. are shown in the parentheses. The contributions to the total \ 2 from each 
experimental data are also shown. 

5 Conclusions 

The AMS-02 collaboration has released the first measurement of the positron fraction 
with unprecedented accuracy. Using the publicly available GALPROP code, we have 
performed a global analysis on the latest data of PAMLEA, Fermi-LAT, and AMS-02, 
in terms of DM annihilation and decay into 2e, 2/i, 2r, 4e, 4/x and At final states. We 
have found that for electron and muon final states, the parameter regions determined 
by AMS-02 is significantly different from that favoured by Fermi-LAT data on the total 
flux of electrons and positrons. The two allowed regions do not overlap at 99.99999% 
C.L.. The conclusion is insensitive to the changes in the normalization and slope of 
the conventional astrophysical background. Consistent fits are obtained only for 2r and 
At final states, which favours TeV scale DM with large cross sections corresponding to 
the boost factor of 0(1000). Such large annihilations are in tension with the current 
measurements of cosmic gamma-rays. Compared with the DM annihilation, the DM 
decay is less favoured by the current data, and so far there is no statistically significant 
indication of charge asymmetric DM decay from the AMS-02 data. 

Note added: As we were finalizing the manuscript, a preprint with similar global 
fitting analysis but different treatment of backgrounds and different focus came out |42j . 
The conclusions in this work are in agreement with theirs. 
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Figure 5: Values of the minimal x 2 from global fits to the data of PAMELA [3], 
Fermi-LAT [4J[8] and AMS-02 [12] as a function of charge asymmetry parameter e for 
DM decaying into 2/i, 4/x, 2r, and 4r final states. The background factors are fixed at 
k = 0.85 and 5 = 0. 
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